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We show that bipolaron formation leads to a d-wave Bose- 
Einstein condensate in cuprates. It is the bipolaron energy 
dispersion rather than a particular pairing interaction which 
is responsible for the d-wave symmetry. The unusual low- 
temperature dependence of the magnetic field penetration 
depth A(T) in cuprates is explained by the localisation of 
bosons in the random potential. Both linear positive and 
negative slopes of A(T) are occur depending on the random 
field profile. 
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The evidence for a (i-like order parameter (changing 
sign when the C11O2 plane is rotated by 7r/2) has been 
reviewed by Annett, Goldenfeld and Legget jjj and more 
recently by Brandow [^| and by several other authors. A 
number of phase-sensitive experiments || provide unam- 
biguous evidence in this direction; furthermore, the low 
temperature magnetic penetration depth |^,^| has been 
found to be linear in many cuprates as expected for a 
d-wave BCS superconductor. However, SIN tunnelling 
studies J@,|7j and some high-precision magnetic measure- 
ments B show the more usual s-like shape of the gap 
function or even reveal an upturn in the temperature 
dependence of the penetration depth below some char- 
acteristic temperature One can reach a compromise 
between conflicting experimental results by mixing s and 
d order parameters or invoking 'anomalous Meissner cur- 
rents' due to surface-induced bound states (often vio- 
lating time reversal symmetry). However, the observa- 
tion of the normal state pseudogap in tunnelling (STM) 
and photoemission spectra (ARPES), non Fermi- liquid 
normal state kinetics and thermodynamics, and unusual 
critical phenomena tell us that many high-T c cuprates 
are not BCS superconductors p0| . In particular, both 
ARPES (ll) and STM @ experiments have shown that 
the maximum energy gap is more than three times larger 
than that expected from the d-wave BCS theory and per- 
sists into the normal state irrespective of doping. The 
gap as well as other major features of STM and ARPES 
have been recently explained with bipolarons Com- 
parison of tunnelling and Andreev gap determinations 
on yittrium, lanthanum and bismuth-based cuprates || 
at various doping levels have unambigously revealed the 
existence of two energy scales as expected for the bipo- 
laronic superconductors |r|. Hence an explanation for 
the d-like order parameter and the anomalous penetra- 



tion depth should be found independent of the BCS gap 
equation. 

While in a BCS superconductor all energy scales and 
symmetries are strictly identical, the symmetry of the 
Bose-Einstein condensate in the bipolaronic supercon- 
ductor should be distinguished from that of the 'inter- 
nal' wave function of a single bipolaron. In this letter 
we show that the Bose-Einstein condensate in cuprates 
is (i-wave owing to the bipolaron energy band structure 
rather than to a particular pairing interaction (see also 
|l2| ) , while the low temperature dependence of the pen- 
etration depth, A(T) is determined by the localisation of 
bipolarons. Both linear positive and negative slopes of 
A(T) occur depending on the random field profile. 

Consideration of perovskite lattice structures shows 
that small inter-site bipolarons are perfectly mobile even 
if the electron-phonon coupling is strong and the bipo- 
laron binding energy is large Jljj ]. Different bipolaron 
configurations can be found with computer simulation 
techniques based on the minimization of the ground state 
energy of an ionic lattice with two holes. The intersite 
pairing of the in-plane oxygen hole with the apex one 
is energetically favorable in the layered perovskite struc- 
tures as established by Catlow et al [ff5|. This apex or 
peroxy-like bipolaron can tunnel from one cell to another 
via a direct single polaron tunnelling from one apex oxy- 
gen to its apex neighbor. The bipolaron band structure 
has been derived by one of us as 

E^ y = tcos(k x , y ) - t'cos(k y , x ). (1) 

Here the in-plane lattice constant is taken as a = 1, t is 
twice the bipolaron hopping integral between p orbitals 
of the same symmetry elongated in the direction of the 
hopping (ppcr) and —t' is twice the hopping integral in 
the perpendicular direction (ppn) . The bipolaron energy 
spectrum in the tight binding approximation consists of 
two bands E x,v formed by the overlap of p x and p y apex 
polaron orbitals, respectively. The energy band minima 
are found at the Brillouin zone boundary, k = (±7r, 0) 
and k = (0, ±7r) rather than at the T point owing to the 
opposite sign of the ppa and ppix hopping integrals. Only 
their relative sign is important, and we choose t, t' > 0. 
Neither band is invariant under crystal symmetry but 
the degenerate doublet is an irreducible representation; 
under a 7r/2 rotation the x band transforms into y and 
vice versa. 

If the bipolaron density is low, the bipolaron Hamil- 
tonian can be mapped onto the charged Bose-gas JUI. 
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Charged bosons are condensed below T c into the states 
of the Brillouin zone with the lowest energy, which are 
k = (±7r, 0) and k = (0, ±7r) for the x and y bipolarons, 
respectively. These four states are degenerate, so the or- 
der parameter ^(m) (the condensate wave function) in 
the real (site) space m = (m x ,m y ) is given by 

*±(m) = N~ 1/2 ^2 b k exp(-ik • m) 

k=(±7r,0),(0,±7r) 

= nl' 2 [cos(7tto x ) ± cos(irm y )] , (2) 

where N is the number of cells in the crystal, 6k is 
the bipolaron (boson) annihilation operator in k space 
(which is a c- number for the condensate), and n c is the 
number of bosons per cell in the condensate. Other com- 
binations of the four degenerate states do not respect 
time-reversal and (or) parity symmetry. The two solu- 
tions, Eq.(||), are physically identical being related by: 
ty + (m x ,m y ) = \P_ (m x ,m y + 1). They have d- wave sym- 
metry changing sign when the CuOi plane is rotated by 
7r/2 around (0,0) for or around (0,1) for (Fig.l). 
The d-wave symmetry is entirely due to the bipolaron 
energy dispersion with four minima at k ^ 0. With the 
energy minimum located at the T point of the Brillouin 
zone the condensate is s-like. 

If the total number of bipolarons in one unit cell is n of 
which jij, are in localised states and no are in delocalised 
states then the number in the condensate n c is 



n c 



n L - n D 



(3) 



and the London penetration depth A cx 1/ ' \fn~ c - Taking 
the delocalised bipolarons to be a free three-dimensional 
gas we have no a Ts. Here we use a simple model we 
have presented previously [[i~7| to calculate the tempera- 
ture dependence of Ul and find that at low temperature 
til varies linearly with temperature. Thus in the limit 
of low temperature we can neglect no and make the ap- 
proximation 



In this limit A(T) — A(0) is small and so 
A(T)-A(0) <xn L (T)-n L (0) 



(4) 



(5) 



i.e. A has the same temperature dependence as til- 

The picture of interacting bosons filling up all the lo- 
calised single-particle states in a random potential and 
Bose-condensing into the first extended state is known in 
the literature Jl^-pT|] . The comprehensive scaling anal- 
ysis of neutral p0| and charged bosons [^lj allows us 
to describe the quantum Bose glass-superfluid transition 
while the thermodynamics of each phase away from the 
transition can be studied with the physically plausible 
models of neutral [^2| and charged bosons jlTj in a ran- 
dom potential. These models are based on a separation 



of localised single-particle states from delocalised states. 
The renormalisation of the single-particle energies by the 
collective mode does not affect the temperature depen- 
dence of any of the thermodynamic functions at low tem- 
peratures |l7| ]. Hence, one assumes that at some temper- 
ature T c bosons are condensed at the mobility edge E c 
so that the chemical potential /i = E Cl and some of the 
bosons are in localised states below the mobility edge. 

For convenience we choose E c = 0. When two or more 
charged bosons are in a single localised state of energy 
E there may be significant Coulomb energy and we take 
this into account as follows. The localisation length £ is 
thought to depend on E via 



l 



-Ef 



(6) 



where v > 0. The Coulomb potential energy of p charged 
bosons confined within a radius £ can be expected to be 



potential energy 



eo£ 



(7) 



where eo is the dielectric constant. Thus the total energy 
of p bosons in a localised state of energy E is taken to be 



w(E) =pE + p{p - 1)k (-Ef 



(8) 



where k > 0. We can thus define an energy scale —E±: 
-Ei = k^. (9) 

From here on we choose our units of energy such that 
Ei = —1. We take the total energy of charged bosons 
in localised states to be the sum of the energies of the 
bosons in the individual potential wells. The partition 
function Z for such a system is then the product of the 
partition functions z(E) for each of the wells, 



where 



p=0 



P0 = -{l + (-E)^} 



and 



(-sy 



-Ei) 



(10) 



(11) 



(12) 



(13) 



The average number of bosons in localised states is 

i-O 



n L = [ dE(p) PL (E) 

J — OO 



(14) 
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where the mean occupancy (p) of a single localised state 
is taken to be 



(?) 



P =oP f 



-a(p-p ) 



p=0 1 



*(p-Po) 



(15) 



and Pl(E) is the one-particle density of localised states 
per unit cell below the mobility edge, which is taken to 
be 



PL\&) = — ei . 
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(16) 



We now focus on the temperature dependence of Ul at 
low temperature (6 -C 1) for the case where the width of 
the impurity tail 7 is large (7 > 1). In the following we 
consider first the case v > 1 and then v < 1. 
If v > 1 we can approximate ul as 



v 



1 



7V L 2(2 - i/) 7 



2^ 



(2-^)7 



In 2 



(17) 



So we expect ul to be close to the total number of wells 
Nl and to increase linearly with temperature. Fig 2a 
compares this analytical formula with accurate numerical 
calculation for the case v — 1.5, 7 = 20. We also note 
that even when 7 < 1, til (6) will still be linear with the 
same slope provided that 9 <C 7. 

If v < 1 we obtain, keeping only the lowest power of 9 
(valid provided 9~ <^ 9) 



_ 1 r(2 - v)^-" 

Nl 2 2 



l-v 



2(2-z/) 7 7 



--In 2 (18) 



Hence in this case Ul decreases linearly with increasing 
temperature (in the low temperature limit). Fig 2b com- 
pares this analytical formula with the numerical calcula- 
tion for the case v = 0.65, 7 = 20. 

Fig 3 shows that the low temperature experimental 
data II on the London penetration depth A of YBCO 
films can be fitted very well by this theory with v < 1 . It 
is more usual to see A increase linearly with temperature 
and this would correspond to v > 1. 

The exponent v depends on the random field profile. 
We believe that v < 1 is more probable for a rapidly 
varying random potential while v > 1 is more likely for a 
slowly varying one. Both v < 1 and v > 1 are observed 
in doped semiconductors. Hence, it is not surprising that 
drastically different low-temperature dependence of the 
London penetration depth is observed in different sam- 
ples of doped cuprates. In the framework of our approach 
A(T) is related to the localisation of carriers at low tem- 
peratures rather than to any energy scale characteristic of 
the condensate. The excitation spectrum of the charged 
Bose-liquid determines, however, the temperature depen- 
dence of A(T) at higher temperatures including an un- 
usual critical behaviour near T c fnOl . 



Many thermodynamic, magnetic and kinetic proper- 
ties of cuprates have been understood in the framework 
of the bipolaron scenario ficj] . We admit, however, that 
one experimental fact is enough to destroy any theory. 
In particular, the single-particle spectral function seen 
by ARPES |ll|] was interpreted by several authors as 
a Fermi liquid feature of the normal state incompatible 
with bipolarons. Most (but not all) of these measure- 
ments indicated a large Fermi surface which one would 
think should evolve with doping as (1 — x) (where x is 
the number of holes introduced by doping) but such an 
evolution is in clear contradiction with kinetic and ther- 
modynamic measurements which show an evolution pro- 
portional to x . Only recently it has been established 
that there is a normal state gap in ARPES and SIN tun- 
nelling, existing well above T c irrespective of the doping 
level 0,UH]. The 'Fermi surface' shows missing seg- 
ments just near the M points j2j| where we expect the 
Bose-Einstein condensation. A plausible explanation is 
that there are two liquids in cuprates, the normal Fermi 
liquid and the charged Bose-liquid (this mixture has been 
theoretically discussed a long time ago [Q). However, 
it is difficult to see how this scenario could explain the 
doping dependence of dc and ac conductivity or of the 
magnetic susceptibility and carrier specific heat. On the 
other hand, the single-particle spectral function of a pure 
bipolaronic system has been recently derived by one of us 
]l2[ . It describes the spectral features of tunnelling and 
photocmission in cuprates, in particular, the tempera- 
ture independent gap and the anomalous gap/T c ratio, 
injection/emission asymmetry both in magnitude and 
shape, zero-bias conductance at zero temperature, the 
spectral shape inside and outside the gap region, tem- 
perature/doping dependence and dip- hump structure of 
the tunnelling conductance and ARPES. The zero-bias 
conductance and any spectral weight at the chemical po- 
tential appear only due to single polaronic states localised 
by the random field. The model is thus compatible with 
the doping evolution of thermodynamic and kinetic prop- 
erties. 

In conclusion we suggest, within the framework of the 
bipolaron theory of cuprates fl(i[| , explanations of the 
d wave symmetry of the ground state and the anoma- 
lous temperature dependence of the London penetra- 
tion depth, compatible with the non Fermi-liquid normal 
state, anomalous critical behaviour, and the ARPES and 
tunnelling data 
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FIG. 2. Dependence of the density of localised bosons til 
on temperature 9 for 7 = 20. (a) v — 1.5, (b) v = 0.65. 
The solid lines correspond to the low temperature predictions 
from equations |l| and [t|, while the dashed lines are derived 
from an accurate numerical calculation. 



FIG. 3. Fit to the London penetration depth obtained by 
Walter et al [8] for a non-irradiated YBCO film. The pa- 
rameter values from the fit were Ei = —74 K, 7 = 20 and 
v = 0.67. 



FIG. 1. D-wave condensate wave 

function, ^-(m) = nl^ 2 [cos(Tvm x ) — cos(Tvm y )] in the real 
(Wannier) space. The order parameter has different signs in 
the shaded cells and is zero in the blank cells. 
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